Measuring the quality factor of a microwave cavity using superconduting qubit devices 
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We propose a method to create superpositions of two macroscopic quantum states of a single-mode microwave 
cavity field interacting with a superconducting charge qubit. The decoherence of such superpositions can be 
determined by measuring either the Wigner function of the cavity field or the charge qubit states. Then the 
quality factor Q of the cavity can be inferred from the decoherence of the superposed states. The proposed 
method is experimentally realizable within current technology even when the Q value is relatively low, and the 
interaction between the qubit and the cavity field is weak. 
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I. INTRODUCTION 

Superconducting (SC) Josephson junctions are considered 
promising qubits for quantum information processing. This 
"artificial atom", with well-defined discrete energy levels, 
provides a platform to test fundamental quantum effects, e.g., 
cavity quantum electrodynamics (QED). The study of the 
cavity QED of a SC qubit, e.g., in Ref. [1], can also open 
new directions for studying the interaction between light and 
solid state quantum devices. These can result in novel con- 
trollable electro-optical quantum devices in the microwave 
regime, such as microwave single-photon generators and de- 
tectors. Cavity QED can allow the transfer of information 
among SC qubits via photons, used as information bus. 

Recently, different information buses using bosonic sys- 
tems, which play a role analogous to a single-mode light field, 
have been proposed to mediate the interaction between the 
SC qubits. These bosonic "information bus" systems can be 
modelled by: nanomechanical resonators (e.g., in Refs. |2|); 
large junctions (e.g., Ref. [3]); current-biased large junctions 
(e.g., Refs. |4]), and LC oscillators (e.g., Refs. |5J|). However, 
the enormous versatility provided by photons should stimu- 
late physicists to pay more attention to SC qubits interacting 
via photons,while embedded inside a QED cavity. 

Several theoretical proposals have analyzed the interaction 
between SC qubits and quantized CI B B S § 113 El or 
classical fields H 1 2i 1 1 3l 1 1 411 . The strong coupling of a single 
photon to a SC charge qubit has been experimentally demon- 
strated 1 15] by using a one-dimensional transmission line res- 
onator 1 16]. But, the QED effect of the SC qubit inside higher- 
dimensional cavities has not been experimentally observed. 
The main roadblocks seem to be: i) whether the cavity qual- 
ity factor Q can still be maintained high enough when the SC 
qubit is placed inside the cavity. Different from atoms, the 
effect of the SC qubit on the Q value of the cavity is not neg- 
ligible due to its complex structure and larger size, ii) The 
higher-dimensional QED cavity has relatively large mode vol- 
ume, making the interaction between the cavity field and the 
qubit not be strong enough for the required quantum opera- 
tions within the decoherence time, iii) The transfer of infor- 



mation among different SC qubits requires the qubit-photon 
interaction to be switched on/off by the external classical flux 
on time scales of the inverse Josephson energy. A higher cav- 
ity Q value, a stronger qubit-photon interaction, and a faster 
switching interaction for the SC qubit QED experiments, seem 
difficult to achieve anytime soon. 

In view of the above problems, it would be desirable to 
explore the possibility to demonstrate a variety of relatively 
simple cavity QED phenomena with a SC qubit. The deter- 
mination of the cavity Q value is a very important first step 
for the experiments on cavity QED with SC qubits. However, 
theoretical calculations of the Q value are not always easy to 
perform because of the complexity of the circuit. Recent ex- 
periments 1 17] on broadband SC detectors showed that the Q 
value of the SC device can reach 2 x 10 6 , which indicates 
that relatively simple experiments using cavity QED with a 
SC qubit are possible. 

In this paper, we propose an experimentally feasible 
method which can be used to demonstrate a simple cavity 
QED effect of the SC qubit. For instance, superpositions of 
two macroscopic quantum states of a single-mode microwave 
cavity field can be created by the interaction between a SC 
charge qubit and the cavity field. At this stage, the injected 
light field is initially a coherent state, which can be easily pre- 
pared. The decoherence of the created superposition states can 
be further determined by measuring either the Wigner func- 
tion of the cavity field or the charge qubit states. Then the 
cavity Q value can be inferred from this decoherence mea- 
surement. Our proposal only needs few operations with a rel- 
atively low Q value. Also, we do not need to assume a very 
fast sweep rate of the external magnetic field for switching 
on/off the qubit-field interaction. Furthermore, the qubit-field 
interaction is not necessarily resonant. 

We begin in Sec. II with a brief overview of the qubit-field 
interaction. In Sec. EH, we discuss how to prepare superposi- 
tions of two different cavity field states under the condition of 
large detuning. In Sec. IV, the cavity Q value is determined 
by the tomographic reconstruction of the cavity field Wigner 
function. In Sec. V, we show an alternative method to deter- 
mine the Q value according to the states of the qubit. Finally, 
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we list our conclusions. 



Hamiltonian Q becomes 



H. THEORETICAL MODEL 

We briefly review a model of a SC charge gubit inside a 
cavity. The Hamiltonian can be written as d 08, 9] 



H = huta? a + E z a z 



(1) 



E,jo x cos 



— ($ C J + 770 + 77* a f ) 
<J>o 



where the first two terms respectively represent the free 
Hamiltonians of the cavity field with frequency uj for the 
photon creation (annihilation) operator (a), and the qubit 
charging energy 



E z = -2E ch (l - 2n g ) , 



(2) 



which depends on the gate charge n g . The single-electron 
charging energy is E c h = e 2 /2(C g + 2Cj) with the capac- 
itors C g and Cj of the gate and the Josephson junction, re- 
spectively. The dimensionless gate charge, n g — C g V g /2e, is 
controlled by the gate voltage V g . Here, a z , a x are the Pauli 
operators, and the charge excited state |e) and ground state \g) 

correspond to the eigenstates | J.) = (?) and | f) ' 



^ 1 J 1 1 ' \0 
of the spin operator a z , respectively. / is an identity oper- 
ator. The third term is the nonlinear qubit-photon interac- 
tion. Ej is the Josephson energy for a single junction. The 
parameter 77 is defined as 77 = J s u(r) • ds with the mode 
function of the cavity field u(r), S is the surface defined by 
the contour of the SQUID. We can decompose the cosine in 
Eq. Q into classical and quantized parts. The quantized parts 
sin[7r(77 a + H.c.) /$o] and cos[7r(?7 a + H.c.) /$o] can be fur- 
ther expanded as a power series in a [a)). To estimate the 
qubit-photon coupling constant, the qubit is assumed to be 
inside a full-wave cavity with the standing-wave form for a 
single-mode magnetic field II 811 



B , 



e Vc 2 



(a — a') cos(fcz). 



(3) 



The polarization of the magnetic field is along the normal 
direction of the surface area of the SQUID, located at an 
antinode of the standing-wave mode. The mode function 
y/ 'hjj I 'sqV c 2 cos(fcz) can be assumed to be independent of 
the integral area because the maximum linear dimension of 
the SQUID, e.g., even for 50 /xm, is much less than 0.1 cm, 
the shortest microwave wavelength of the cavity field. Then, 
in the microwave regime, the estimated range of values for 
7r?7/$o is: 8.55 x 10~ 6 < vr77/$ < 1.9 x 10~ 3 , for a fixed 
area of the SQUID, e.g., 50 /im x50 /jm. If the light field is 
not so strong (e.g., the average number of photons inside the 
cavity N = (aJa) < 100), then we can only keep the first 
order of ir-q/Qo and safely neglect all higher orders. Thus, the 



H = huja)a + E z a z — Eja x cos(— — ) 

<t>o 



TIE.] , 7T$ 

sm 

$0 *o 



sm( — ) (7700+ +77 a 1 cr_) 



(4) 



It is clear that the qubit-photon interaction can be controlled 
by the classical flux $ c , after neglecting higher-orders in 

7T?7/$o- 



III. PREPARATION OF MACROSCOPIC 
SUPERPOSITION STATES 



The qubit-photon system can be initialized by adjusting the 
gate voltage V g and the external flux <1> C such that n g = 1/2 
and $ c = 0, then the dynamics of the qubit-field is governed 
by the Hamiltonian 



E T a 



(5) 



Now there is no interaction between the cavity field and the 
qubit; thus, the cavity field and the qubit evolve according to 
Eq. 0. We assume that the qubit-photon system works at 
low temperatures T ( e.g., T — 30 mK in Ref. 0), then the 
mean number of thermal photons (nth) in the cavity can be 
negligible in the microwave regime |8], and the cavity is ap- 
proximately considered in the zero temperature environment. 
The initial state of the cavity field is prepared by injecting a 
single-mode coherent light 



exp 



(6) 



into the cavity. Here, without loss of generality, a is assumed 
to be a real number, and a\a) — a\a). The qubit is assumed 
to be initially in the ground state \g). After a time interval 
Ti = fm/AEj, the qubit ground state \g) is transformed as 
1.9) ~~ * (Iff) + *l e )) / v2; then, the qubit-photon state evolves 
into 



1 



|^(ri)> = -^(| 5 )+i|e) 



(7) 



Here we have neglected the free evolution phase factor e~ iuJTl 
in a. 

Now, we assume that the gate voltage and the magnetic flux 
are switched to n g ^ 1/2 (this value of n g will be specified 
later) and $ c = $o/2, respectively. Then the qubit-photon 
interaction appears and the effective Hamiltonian governing 
the dynamic evolution of the qubit-photon can be written as 
(see Appendix A) 



Ho = huj-a/a 



1 



A 



;i + 2at a ) | e )( e |, (8) 



with w_ = uj - |g| 2 /A and g = (nr]Ej)/(h$o). The detun- 
ing A = — uj > between the qubit transition frequency 
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n = —4:E c h(l — 2n g )/h and the cavity field frequency to is 
assumed to satisfy the large detuning condition 



< 1. 



The unitary evolution operator corresponding to Eq. 
be written as 



n 



i uj-d'a + —<7 Z t 



U(t) = exp 

x exp [— itF{a) a)\e){e\\ , 
here, the operator F(a^a) is expressed as 



(9) 



(10) 



(a) 



IT 




FIG. 1: (Color online) Wigner functions W+(x,p) of Eq. i 14i and 
W+d(x,p) of Eq. 1 16t for the cavity field without and with the en- 
ergy dissipation are shown in (a) and (b), for the input state |a) with 
n — 16. Here the Wigner functions W+(x,p) and W+d{x,p) are 
normalized to n N+ . 



F(a f a) = J^-(l + 2a + o). (11) 
With an evolution time r 2 , the state (0 evolves into 

fo&fa)) = ±= \\ g )\{3) +iex P (^)|e)|/3')] (12) 

where a global phase exp(— if2r 2 /2) has been neglected, 9 = 
(fi — |g| 2 /A)r 2 , /3 = aexp[— iw_r 2 ], and (3' — [3 exp(— i<j>), 
with (f> = 2|g| 2 r 2 /A. Equation d!2t shows that aphase shift 
is generated for the coherent state \0) of the cavity field when 
the qubit is in the excited state |e), but the qubit ground state 
\g) does not induce an extra phase for the coherent state \(3). 

The gate voltage and the magnetic field are now adjusted 
such that the conditions n g = 1/2 and $ c = are satisfied; 
then the qubit-photon interaction is switched off. Now let the 
system evolve a time r' = t% — fm/AEj, then Eq. dl 2i be- 
comes 



|V>(t 2 )) = \\g) 
■ 1, \ 

+ *- 2 \e) 



hexpW)]?)], (13) 



where a free phase factor e~ WTl in the cavity field states \(3) 
and |/3e~ 4 ^) has been neglected. 

The superpositions of two distinguished coherent states can 
be conditionally generated by measuring the charge states of 
the qubit as, 



\fti 



N ± l [ 



(14) 



where the + (— ) correspond to the measurement results |e) 
(\g)), and the normalized constants N± are determined by 



Nl = 2 ± 2 cos 9' exp 



-2|a| 2 sin 2 



(15) 



where 9' — \a\ 2 smtfi— 9, and the relation \/3\ 2 = \a\ 2 is used. 

Due to $ c = 0, after the superpositions in Ea. (ll4> are cre- 
ated, the dynamic evolution of the cavity field is only affected 
by its dissipation, characterized by the decay rate 7, which 
can be expressed by virtue of the cavity quality factor Q as 
Q = lli/j. Now let the cavity field described by the states 



) or evolve a time T3; then the reduced density matri- 
ces of the superpositions can be described by 



P±{tz) = -^ i {\Pu)(Pu\ + \13 i u)(P'u\ 

iV ± 

±C\(3'u)(f3u\±C*\f3u)((3'u\}, 



(16) 



where 



and 



C = exp(i6»)exp{|a| 2 (l-e-^)(u 2 - 1)} (17) 



u{t 3 ) = exp (-7^3) 



exp 



WT3\ 

2Qj 



(18) 



It is clearly shown that the mixed state in Eq. fl!6l > is strongly 
affected by the Q value. Equation dl6l is derived for zero 
temperature since thermal photons are negligible at low- 
temperature. Equations dl6!18l > show that the information of 
the cavity quality factor Q can be encoded in a reduced den- 
sity matrix of the cavity field. The Q value can be determined 
using two different methods, after encoding its information in 
Eas. (ll6ll8l >. Below, we will discuss these two approaches. 



IV. MEASURING Q BY PHOTON STATE TOMOGRAPHY 

The state p of the optical field is generally visualized when 
it is represented by a Wigner function |29] in the position x 
and momentum p space, which is written as 



1 f°° 

W(x,p) = - / (x- x'\p\x + x')e l2px 'dx'. 



(19) 



The Wigner function W(x,p) can be experimentally mea- 
sured by state tomographic techniques 1 29] . For any two co- 
herent states, \a) and |/3), the Wigner function W(x,p) can 
be represented as 



W a ,p{x,p) 



(x- x'\a)((3\x + x')e l2px 'dx' 
1 



■ exp 



-(H 2 + |/?| 2 -2a/r) 



exp {-(x - qi) 2 -(p+ iq 2 ) 2 } 



(20) 
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with q-L = (a + (3*)/V2 and q 2 = (a - j3*)/\f2. The 
Wigner functions W±{x,p) and W±d(x,p) for the states 
(114-i and ( 1161 were calculated (see Appendix B) by using 
Eq. J20I . Comparing the tomographic ally measured results 
for the states (114-i and d!6i . the Q factor of the cavity can be 
finally determined, as explained below by using an example. 

We further numerically calculate the Wigner functions 
W±(x,p) and W±o(x,p) of the states dl4t and d!6l > from 
the SC qubit parameters and given operation durations. Us- 
ing current values for experimental data, the basic physical 
parameters can be specified as follows. We assume that the 
SC Cooper-pair box is made from aluminum, with a BCS en- 
ergy gap of ~ 2.4K (about 50 GHz) 12011 . the charge energy 
_E c h and the Josephson energy Ej are AE^/ h = 149 GHz and 
2Ej/h — 13.0 GHz, respectively |20j. The frequency of the 
cavity field is taken as 40 GHz, corresponding to a wavelength 
~ 0.75 cm. The above numbers show that the SC energy gap 
is the largest energy, so the quasi-particle excitation on the is- 
land can be well suppressed at low temperatures, e.g., 20 mK . 
The SQUID area is assumed to be about 50 /im x 50 /im, then 
the absolute value \g\ of the qubit-photon coupling constant is 
about \g\ = 4 x 10 6 rad s _1 . 

Let us now prepare entangled qubit-photon states as in 
Eq. (II 31 . Any gate charge value n g in Eq. (0, in which the 
large detuning condition in Eq. (|9} is satisfied, can be chosen 
to realize our proposal. For concreteness, we give an exam- 
ple. The gate voltage is adjusted such that the gate char ge i s 
n g ~ 0.634233, which can be experimentally achieved 1 20], 
then the detuning A = fl - w « 9.0 x 10 6 rad s _1 . Thus, 
fl is about 40 GHz plus 1.4 MHz, and \g\ 2 /A ~ 0.27 MHz. 
We can also find that A/lfll ~ 2.3, so a large-detuning con- 
dition can be used ll2ll 12211 . For a given Josephson energy, 
2Ej/h = 13.0 GHz, the operation times t x = 4.8 x 10~ 12 
s, required to prepare a superposition of |e) and \g) with 
equal probabilities, is much less than the qubit relaxation time 
T\ = 1.3 /is and dephasing time T2 = 5 ns. We can choose 
the duration T2 for a given input coherent state |a) with the 
condition, that the distance \/3 — f3'\ between two coherent 
states |/3) and |/3') satisfies 

\f3-P'\ = 2Hsin(^ > 1. (21) 

So the lower bound of the duration T2 can be given as 

T2 = A. arcsin (J_) , (22 ) 

when < cj) < ir. Equation d22i shows that a shorter r 2 can 
be obtained for a higher intensity \a\ with fixed detuning A 
and coupling constant g. 

As an example, we plot the Wigner function of the super- 
position |/3+) in Fig.^a) for an input coherent light |a) with 
the mean photon number n = \a\ 2 = 16. We choose a simple 
case (j> = tt, corresponding to the operation time T2 ~ 0.93 /is, 
which is less than the qubit lifetime T\ and the cavity field 
lifetime T p ^ « 2 /is for a bad cavity with Q — 5 x 10 5 . In 
such a case, (3' = —(3 and the phase 8 is about 0.996 [mod 2ir] 
rad. Other parameters used in Fig. [2 are given above. If we 



set the evolution time T3 = 0.1 /is, then the Wigner func- 
tion of Eq. ( I16> for the above cavity quality factor is shown 
in Fig. 0b). The central structure in Fig.Qa) represents the 
coherence of the quantum state. In Fig.^b), we find that the 
height of the Wigner function W+d(x,p), especially for the 
central structure, is reduced by the environment. Comparing 
Fig. Qa) and Fig. \l[b), it is found that the coherence of the 
superposed states is suppressed by the environment, and the 
decoherence of superpositions is tied to the energy dissipation 
of the cavity field. Then, the Q value can in principle be es- 
timated by measuring the Wigner functions of Eqs. (116b and 
(I14i . and comparing these two kinds of results. 



V. DETERMINING Q BY READOUT OF CHARGE STATES 

The determination of the Q value by measuring the Wigner 
function needs optical instruments. In solid state experiments, 
the charge states are typically measured. Instead of using op- 
tical instruments, it would be desirable to obtain the Q value 
by measuring charge states. This will be our goal here. The 
process to achieve this can be described as follows. 

i) According to the measurements on the charge qubit states 
in Eq.(I13>. the qubit-photon states are projected to \g) Cg) |/3_) 
or |e) ® |/3+), respectively. After the evolution time 7-3, a tt/2 
quantum operation is performed on the qubit with the duration 
T = hn/AEj. Then, the qubit ground state \g), or excited 
state |e), is transformed into the superposition (|<?)+i|e))/ V2, 
or (i\g) + |e))/v2, and the photon states \/3±) evolve into 
mixed states after the evolution time r = T3 + T, and the 
photon-qubit states can be expressed as 

PQ+f = \{\g) ± i|e»«ff| T i(e\) ® p±(t) , (23) 

with subscripts Q and F denoting the qubit and cavity field, 
respectively. The reduced density matrices p±(r) take the 
same form as in Eq. (I16> with r replacing T3. 

ii) After the above procedure, the qubit-photon interaction 
is switched on by applying the external magnetic flux $ e = 
$o/2. By using Eq. J10I . Eq. (I23> evolves into 

IPA+F = \9){9\®U 1 {T i ) P ±{T)Ut{T i ) 

+ \e)(e\®U 2 (T 4 ) P± (T)ul(T 4 ) (24) 
T zexp(-ift_T 4 )|5)(e| ® U 1 (t 4 )p±(t)uI(t 4 ) 
± iexp(-MO_T 4 )|e)(#| <x> U 2 (t 4 )p±(t)uI(t 4 ) 

with fl- = CI — |g| 2 /A, and a shorter evolution time t 4 . For 
example, t 4 is less than the lifetime T\ of the qubit at least. 
The time evolution operators Ui(t 4 ) and U2(t 4 ) in Eq. (1241 
are 

Ui(t 4 ) = exp ^074] , (25a) 

^2(74) = cxp [— iuj + a^a T4] . (25b) 

with ui± = ijj ± |g| 2 /A. After this qubit-photon interaction, 
the information of the Q value is encoded. 
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iii) The qubit-photon coupling is switched off and a 7r/2 
rotation is made on the qubit. If the state of the cavity field is 
prepared to |/?_) of Eq. (114-1 in the first step, then the qubit is 
in the ground state \g). After measuring the qubit states, the 
photon states are projected to 



Pe/g 



\{A±B) 



(26a) 



where the sign " + " corresponds to the excited state |e) mea- 
surement, but " — " corresponds to the ground state \g) mea- 
surement. The operators A and B are 

2 

A = Y,Ui(T 4 ) P -(T)U}(T4), (26b) 
i=l 

B = 2Re[exp(-^_r4)C/i(r4)p_(T)C/ 2 t (r 4 )]. (26c) 

After tracing out the cavity field state, the probabilities cor- 
responding to measuring charge states |e) and \g) are 



Pe/ g (r) 



Tr F (p e/g 
1 



{l±Re(Tr F [exp(-^)p_(r)])} (27) 



with ip = — 2\g\ 2 a' f a/A)T4. Then the measurement prob- 
abilities are related to the Q values. Substituting P-(t) into 
Eq.(l27L we can obtain 



Re{Tr[exp(-iy>)p_(r)]} 
2 



(28a) 



= exp [— 2a(r) sin 2 </>'] cos [J7_r 4 — a(r) sin(2</>')] 
--^j cos [6L - |a| 2 sin</> + 9 - fi_r 4 ] exp (+G_ - Y) 
cos [Q+ + l a l 2 sin0 - 6» - fi_r 4 ] exp (-G+ - T) 
with the parameters 





\9\\ 




(28b) 


r 


= 2|a| 2 sin 2 ( 


2 j ' 


(28c) 


a(r) 
G± 

o± 


= Mr)| 2 , 
= 2a(r)sin(/ 
= 2a(r) cos( 


/ sin(</»±0'): 
(j>±</)')sm <p'. 


(28d) 
(28e) 
(28f) 



From Eq. d28i . we find that </>' should satisfy the condition 
4>' 7^ rm for = 7r, in order to describe the dissipation effect; 
here n is an integer. Generally speaking, if one of the func- 
tions G+, G_, 9 + , 9-, sin <fi', or sm(2<fi') is nonzero, then this 
is enough to encode the Q value, which can be obtained from 
Eq. (I28> . together with Eq. Jl 81 . using r instead of T3. 

However, if the superposition of the cavity fields is prepared 
to the state in the first step, then the ground and excited 
state measurements make the cavity field collapse to state 




FIG. 2: (Color online) The probability P g (r) to measure the qubit 
ground state \g) as a function of the evolution time r. This Pg(r) 
is shown for several values of the quality factor Q and for different 
intensities of the input coherent state |a). 



where A' and B' have the same forms as Eqs. J26bl and ( I26ct . 
just with the replacement of p_ (t) by p + (r). 

The probabilities Pg(r) and P' e (r) to measure the qubit 
states \g) and |e) corresponding to the prepared state \(3 + ) of 
Eq. dl4> after a dissipation interval t, can also be obtained as 



PL 



1 



g/e (r) = - {1 ± Re(Tr F [exp(-^)p + (r)])} (30) 



where Re{Tr F [exp(— iip)p + (r))} can be obtained by replac- 
ing N_ with N + , and replacing the sign " — " before the sec- 
ond and third terms with the sign " + " in Eq. ( 128b ) 

To determine the Q values by probing the charge states, the 
measurement should be made for two times, the first measure- 
ment is for the preparation of the superpositions of the cavity 
field. After the first measurement, we make a suitable qubit 
rotation, and then make the qubit interact with the cavity field 
for a duration r 4 . Finally, the second measurement is made 
and the Q information is encoded in the measured probabili- 
ties. The different evolution times r correspond to the differ- 
ent measuring probabilities for given r 4 and other parameters 
\a\, A, and so on. For example, the probabilities P e / ff (r) for 
several special cases are discussed as follows when the pre- 
pared state is |/3-). If we assume that the qubit rotations and 
qubit-photon dispersive interaction are made without energy 
dissipation of the cavity field, e.g., r = 0, then the measuring 
probabilities P e / S (r = 0) only encode the information of the 
cavity field but do not include the quality factor Q. If the co- 
herence of the states \/3±) nearly disappears after time r, then 
the state |/3_) becomes a classical statistical mixture 



P-(j) = [\Mt))(Mt)\ + \(3'u(t))((3'u(i 



(31) 



(32) 



Pg/e 



(29) 



The probabilities P e / g (r) are then reduced to 

J_ exp [-2a(-r) sin 2 <jf\ 

r e/g(T) - N 2 ± N 2 

x cos [Q_r 4 — a(r) sin(20')] , 



which tends to 1/2 for \a\ 2 > 1. If r > I/7 = t ph of 
the single-photon state lifetime, then the photons of the states 
\{3±) are completely dissipated into the environment. In this 
case, the cavity quality factor Q cannot be encoded in the 
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probabilities P e / g (r) even with some qubit and qubit-photon 
states operations. 

As an example, let us consider how P g (r) varies with the 
evolution time r with the cavity field dissipation. We assume 
that the evolution time r 4 = (7r/2)(A/|g| 2 ), that is, <fi' = n/2. 
Then, the r-dependent probabilities P g (r) for the initially pre- 
pared state |/3_) are given in Fig.|2ja) with the same param- 
eters as in Fig.^ except with different cavity quality factors 
Q. In order to see how the probability P g {r) changes with the 
intensity \a\ 2 of the input coherent state \a), we plot P g (r) in 
Fig- 0b) with the same parameters as Fig.|2fa) except chang- 
ing the intensity to \a\ 2 = 4 from \a\ 2 = 16. Figure [2] shows 
that both the higher quality factor Q and weaker intensity \a\ 2 
of the input cavity field correspond to a larger probability P g 
of the ground state for the fixed evolution time t. For fixed 
Q and t, the weaker intensity \a\ 2 corresponds to a higher 
measuring probability. We plot P g (r) in Fig. [2] considering 
the simple case <f> — it. However, if we consider another <f>, 
then | a | 2 should be chosen such that it satisfies the condition 
2\a\ sm(cf>/2) > 1. In conclusion, the quality factor Q can be 
determined from the probabilities P e / g (r) of measuring the 
qubit states with a finite cavity field evolution time r. 



VI. DISCUSSIONS AND CONCLUSIONS 



Although we did not give a detailed description of another 
resonance-based approach, it should be pointed out that the Q 
values can also be determined by virtue of the resonant qubit- 
photon interaction. For example, if the superpositions |8] of 
the vacuum and the single photon state are experimentally pre- 
pared, then we can follow the same steps as in Sec. Ill and IV 
to obtain the Q value. This method [24] has been applied to 
micromasers, where the qubits are two-level atoms. However, 
the coherent states and non-resonant qubit-photon interaction 
should be easier to do experimentally than the approach using 
single-photon states and resonant qubit-photon interaction. 

Our proposal can also be generalized to the models used 
in Refs. I15l 1 1611 . which are experimentally accessible. We 
hope that our proposal can open new doors to experimentally 
test the Q value and motivate further experiments on cavity 
quantum electrodynamics with SC qubits. 
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We discussed how to measure the cavity quality factor Q 
by using the interaction between a single-mode microwave 
cavity field and a controllable superconducting charge qubit. 
Two methods are proposed. One measures the Wigner func- 
tion of the state d!6t by using a standard optical method \ 29]. 
Another approach measures the qubit states. Using this last 
method, the information of the Q value can be encoded into 
the reduced density matrix of the cavity field, and at the same 
time the qubit makes a ir/2 rotation. Thus, with a suitable 
qubit-photon interaction time, information on the Q value is 
then transferred to the qubit-photon states. Finally, after an- 
other tt/2 rotation, the charge qubit states are measured, and 
the Q value can be obtained, as shown in Fig. [2J Eqs. fl28i 
and Eq. dl8l l. However, it should be noticed that it is easy to 
measure charge states than to measure photon states in super- 
conducting circuits. 

Our proposal shows that a cavity QED experiment with a 
SC qubit can be performed even for a relatively low Q values, 
e.g., Q ~ 10 6 . Initially, a coherent state is injected into the 
cavity, which is relatively easy to do experimentally. Although 
all rotations of the qubit are chosen as w/2 to demonstrate our 
proposal, other rotations can also be used to achieve our goal. 

To simplify these studies and without loss of generality, we 
have assumed two components \(3) and \/3 exp(i(f>)) for su- 
perpositions with (f) — n phase difference in our numerical 
demonstrations. Of course, other superpositions can also be 
used to realize our purpose. The only condition to satisfy is 
that the distance between the two states \(3) and |/3exp(i0)) 
should be larger than one. In order to obtain a numerical es- 
timate for the detuning, we specify a value of the gate charge 
number n g . However, any gate charge that satisfies the large- 
detuning condition can be chosen to realize our proposal. 



APPENDIX A: EFFECTIVE HAMILTONIAN WITH 
LARGER DETUNING 

The Hamiltonian H = Hq + Hi of the two-level atom in- 
teracting with a single-mode cavity field can be written as 

Hq = —hVl<j z + hua^a, (Ala) 

Hi = h(ga,i(T- + g*a<r+) (Alb) 

with a complex number g. Let us assume A = SI — u> > and 
gj A <C 1. The eigenstates and corresponding eigenvalues of 
the free Hamiltonian Hq are 

|e) ® \ri) nhu + -Ml, (A2a) 

\g) (g) | to) => mfuv - (A2b) 

In the interaction picture, any state can be written as 

\^(t)) = U(t,t )m )) (A3) 

with 

1 I* 

Ufa to) = 1 + / ffint(*i)dti (A4) 

+ (J-J £ £ H int {ti)H int {t 2 )dtidt 2 + --- 

here H int = C/ t (t)i/ 1 C/ (i) with U (t) = exp{-iH t/h}. In 
the basis {\Ei}} = {|e) ® |n), \g) <g> |m)}, Eq. dA4t can be 
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expressed as 

[/(Mo) = 1 + 
1 



with 



ih J t 



£)|£l}<£j|iW*l)|£m>CEm|d*l 



l.m 



After neglecting the fast-oscillating factor and keeping the 
first order term in g /A, U (t, to) 



U(t,t ) = U(t,0) = U(t) 



(A5) 



wl-i— y d*i[(n + l)|e,n)(e,n|-n|fli,n)(5,n|], 

where we assume to = 0. Finally, we obtain the effective 
Hamiltonian in the interaction picture as 



H, 



off 



■{\e){e\aa^ - \g){g\<Ja) 



(A6) 



Returning Eq. ( IA6> to the Schrodinger picture, Eq. dlOi is ob- 
tained. This method can be generalized to obtain the effective 
Hamiltonian of the model with many two-level system inter- 
acting with a common single-mode field. Equation (II 01 can 
also be obtained by using the Frohlich-Nakajima transforma- 
tion |0 nielli. 



APPENDIX B: WIGNER FUNCTIONS OF SUPERPOSITION 
AND MIXED STATES 

For completeness, we explicitly write the Wigner functions 
W±(x,p) of the superposition states in Eq. (II 41 as follows: 



W±(x,p) 
1 



*N* 1 6XP 



+ exp 
±2Re 



(x-V2Re/3) -(p-V21mj3 
Pexp(-(.T- Pl ) 2 -(p + ip 2 ) 2 )]} , (Bl) 



P = 
Pi = 

P2 = 



expH60ex P Ha| 2 (l-e^) 



(B2) 
(B3) 

(B4) 



The Wigner functions W±d(%,p) of the mixed states in 
Eq. ( 1161 with dissipation can be written as 



W±d(x,p) 
1 



irNl 



+ exp 
±2Re 



(x - uV2Ref3^j -(p-uV21mp 
- (x - u%/2Re/3') 2 - (p - uV2Im(3'y 
Pexp (- (x- upif - (p + iup 2 ) 2 ^J |. (B5) 
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